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We study the Kelvin mode excitations on a vortex line in a three-dimensional trapped Bose-Einstein conden- 
sate at finite temperature. Our stochastic Gross-Pitaevskii simulations show that the activation of these modes 
can be suppressed by tightening the confinement along the direction of the vortex line, leading to a strong 
suppression in the vortex decay rate as the system enters a regime of two-dimensional vortex dynamics. As 
the system approaches the condensation transition temperature we find that the vortex decay rate is strongly 
sensitive to dimensionality and temperature, observing a large enhancement for quasi-two-dimensional traps. 
Three-dimensional simulations of the recent vortex dipole decay experiment of Neely et ah [Phys. Rev. Lett. 
104, 160401 (2010)] confirm two-dimensional vortex dynamics, and predict a dipole lifetime consistent with 
experimental observations and suppression of Kelvon-induced vortex decay in highly oblate condensates. 



I. INTRODUCTION 

Quantum vortices in Bose-Einstein condensates (BECs) are 
topologically stable excitations of the superfluid matter wave 
field that carry angular momentum ^ . Unless the system is 
rotatingsufficiently fast, vortices are thermodynamically un- 
stable 1121. At finite temperatures the thermal cloud can pro- 
vide a source of dissipation and noise, and hence allow vortex 
decay. However, there is currently little understanding of the 
decay mechanism and its timescale. The dynamics of vor- 
tices in such systems is therefore of interest as a robust test of 
non-equilibrium field theory. A number of different methods 
have been developed to describe the dynamics of vortices in 
finite-temperature systems. Single vortex decay has been sim- 
ulated using the Zaremba-Nikuni-Griffin formalism [jsl], the 
classical field method |4], and within a Gross-Pitaevskii the- 
ory with phenomenological damping The nature of the 
thermal instability has also been studied within a variational 
approach [6] to the stochastic Gross-Pitaevskii equation [7]. 

Classical vortices can support long-wavelength helical trav- 
eling waves known as Kelvin waves | 8]. Superfluid Kelvin 
waves, or Kelvons |9], have long be studied in Helium ifToll . 
and more recently have become a topic of interest in BECs. 
The Kelvon excitation mechanism has been theoretically in- 
vestigated ifTTI] . along with the microscopic dynamics of co- 
herently excited Kelvin waves with results consistent 
with the experimental evidence for Kelvons |[T3|] . Similar 
bending waves have also been shown to play an important role 
in the instability of vortex rings in trapped BECs [Il4ll . 

While the coherent excitation of Kelvin modes in a BEC 
is well understood, their role in the finite-temperature system 
is largely unexplored. Here we show that the thermal activa- 
tion of Kelvin waves is a dominant factor in the decay of a 
vortex for a finite-temperature three-dimensional (3D) BEC. 
The underlying reason is that Kelvin waves cause the vortex 
to wobble and emit acoustic radiation, thus allowing the vor- 
tex to more effectively dissipate energy, and hence decay by 
moving out to the boundary of the BEC. As Kelvin waves are 
fundamentally 3D excitations of the vortex line, they may be 
suppressed by flattening the system (see Fig.[T]). 

Indeed, beyond a certain oblateness we find that the Kelvon 
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FIG. 1. Schematic of vortex (solid line) bending from the z-axis 
(dashed line) in a trapped BEC for (a) a spherical trap (in terms of 
trap frequencies defined in Eq.|3] K = uJz/uJr ^ 1), and (b) a highly 
oblate trap (A > 1). 



mode contribution is essentially frozen out and the vortex de- 
cay rate reduces to a geometrically-invariant value, character- 
izing a regime of two-dimensional vortex (2DV) dynamics. 
We also find that for tighter traps and higher temperatures the 
matter wave itself crosses over to being quasi-2D and we find 
a dynamical signature of the phase-fluctuating condensate: an 
anomalous increase in the rate of vortex decay. 

As a further application of these concepts, we model a 
recent experiment [Hsll studying vortex dipoles in the 2DV 
regime, confirming the 2DV behavior in finite temperature 
dynamics. The decay characteristics are found to be quali- 
tatively very different from 3D dipole decay which exhibits 
vortex ring and loop formation, accelerating the damping pro- 
cess. 

Our theoretical framework for studying this system is the 
stochastic projected Gross-Pitaevskii equation (SPGPE) llT6ll 
which is a grand-canonical c-field method jvh (also see lIlSl 
19] for related theories). This formalism has been used to 
study vortex lattice formation from a rotating thermal cloud 
||2Q|i , a quantitative model of spontaneous vortex formation in 
experiments JIO, lUl], and the decay of a single vortex ||22|] 
(also see 1I2I ). 
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FIG. 2. Density isosurfaces of the C-field after t = 0.61 s of SPGPE evolution, for a range of trap geometries. Red: the full density isosurface. 
Blue: isosurface restricted to the region of the vortex core. The system temperature is T = 0.78Tc . 



II. SYSTEM AND EQUATION OF MOTION 

A. Physical system 

We begin by describing the system of bosonic atoms in the 
cold-collision regime with the second-quantized many body 
Hamiltonian 

+ 1 J d?x ^(x)'^^(x)'^^(x)^(x) (1) 



where the single-particle Hamiltonian is 



2m 



the cylindrically symmetric harmonic trap is 



m 



(2) 



(3) 



and g = ATrTi^a/m, with m the atomic mass and a the s-wave 
scattering length. The field operators obey Bose commutation 
relations 



[^(X),^t(y)]=^(x-y). 



(4) 



We define A = cj^/cj^ as the degree of trap oblateness, with 
uj^ = A~^/^a), enforcing a constant geometric mean fre- 
quency Cd^ = uj^uoz for arbitrary oblateness. This choice en- 
sures the ideal gas transition temperature is invariant under 
scaling of A at fixed total atom number N. 



B. Stochastic PGPE Formalism 

In this work we use the SPGPE formalism to describe the 
process of vortex decay in a finite temperature BEC. The ap- 
proach we follow in choosing parameters and preparing initial 
states for the SPGPE is extensively detailed in Ref. | 22], hence 
we only briefly summarize it here. 

The SPGPE is a c-field method where the system is de- 
composed in terms of the single-particle modes of the system 



satisfying i^sp0n(x) 



(x), where n represents all quan- 



divided into the coherent region (C) consisting of all states 
with Cn < Ccut, and the incoherent region (I) containing the 
remaining high energy modes [17]. This division is made us- 
ing a projection operator, Vq which defines the C-region field 
operator as 



7/)(x) = Vc^(yi) 



C(y)*(y)- (5) 



Equivalently, we have the expansion 



V^(x) = E ^^n0n(x), 



(6) 



where [a^, aj^] = 5n,m- The I-region acts as a thermal reser- 
voir for the C-region, and is assumed to be in thermal equi- 
librium at a temperature T and chemical potential ji. The 
C-region dynamics are described by a stochastic differential 
equation for the c-field corresponding to ©I 



(7) 



where an are ordinary c-numbers. The equation of motion is 
known as the simple-growth SPGPE ||2Q|] 



{^-Lc)^cdt^dw], (8) 



where Lq is the Hamiltonian evolution operator for the C re- 
gion LcV^c = (^sp + ^IV^cP) V^c. and the complex Gaussian 
noise, dW, satisfies (6iW*(x, t)^^W(x^ t)) = 2-f6{^-id)dt, 
(diy (x, t)dW{-x',t)) = 0. For a fufl derivation of the SPGPE 
we refer the reader to [ 16] . 

The first term of the SPGPE ^ describes Gross-Pitaevskii 
evolution over the C region, while the second and third terms 
account for scattering of two I-region atoms, resulting in 
growth of the C-region, and the corresponding time-reversed 
process. The growth rate, 7, can be calculated in near- 
equilibrium systems and is independent of position over 
a large region of the trap with the value 
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tum numbers required to specify the modes. The system is where p = l/Zc^T and <l> is the Lerch transcendent. 
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FIG. 3. Power spectra for a range of geometries as a function of 
kz/k, where k = n /W corresponds to the largest wavelength re- 
solvable for the chosen window W — Rtf 1 2 • The spectra are 
calculated after t = 0.70 s, and for T = 0.78T°. 



To obtain initial states for evolution according to (|8]), we 
follow the procedure used in [I22I1 . which we briefly summa- 
rize. We generate (non-vortex) finite-temperature equilibrium 
states by evolving the SPGPE with an appropriate choice of T, 
/i, and CcutiHl- The values of T, /i, and ecut then determine 
7 [Eq. (|9])] and hence all SPGPE parameters are obtained in a 
physically consistent manner prior to simulation. To create the 
vortex initial state we imprint a phase V^c(x) ?/^c(x)e*®*^^^ 
onto the equilibrium state, where 6(x) is the azimuthal angle 
in the radial plane [25]. This creates a straight singly-charged 
vortex line along the z axis of the trap. This procedure gener- 
ates a non-equilibrium state because the reservoir parameters 
in dSj) are for a non-rotating reservoir. 



III. SINGLE VORTEX DECAY 

In our simulations we use initial vortex states prepared for a 
system containing 5x10^ ^^Rb atoms in a trap of (fixed) geo- 
metric mean trap frequency Gj = 27r x 19.7s~^ [for which 
= 69.7nK], and a range of geometries from spherical 
to highly oblate. For each parameter set we calculate vor- 
tex decay properties by evolving 10-30 trajectories of the 
SPGPE iImI. 

Fig. [2] shows density isosurfaces of the vortex in a va- 
riety of condensate geometries after t = 0.61 s of evolu- 
tion. The phase imprinting technique leads to a straight vortex 
line along the z-axis, but under evolution the vortex rapidly 
thermalizes, typically within 0.1 s (see discussion below). 
We observe a significant amount of vortex bending in the 
spherical geometry and as the condensate becomes more two- 
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FIG. 4. (a) Total power as a function of trap oblateness A. (b) Ef- 
fective vortex decay rate (1/t) as a function of A. The disconnected 
points for A = 20 break from the trend of vortex decay rate indepen- 
dent of A, as discussed in section UlTcl 



dimensional, the magnitude of bending lessens and the vortex 
behaves essentially as a straight line. Although vortex bending 
on a short length scale is evident in spherical geometries, the 
vortex stays approximately aligned with the z-axis throughout 
its evolution. 



A. Kelvon power spectrum 

To obtain the vortex line coordinates during evolution we 
locate the radial position of the phase singularity {xy^y^) in 
each plane (i.e. z value) and for all times t. Fluctuations in 
the phase at the condensate boundary restrict the spatial range 
over which the vortex can be detected. This reduces the range 
of z-values for which we can obtain the vortex-line coordi- 
nates, and, as time evolves and the vortex precesses toward 
the radial condensate boundary, this range of z-values further 
reduces. For this reason we limit vortex detection to the range 
z = ±W, where W = ^^tF' ^i^^ ^tf Thomas-Fermi 
radius along the z direction. We have found this to be the 
largest range which gives a well defined and unique curl sig- 
nal over the vortex lifetime for our simulations. 

We define 2i(z, t) = [xy{z^t)—x{t)]-\-i[yy{z^t)—y{t)] to he 
the (complex) vortex line displacement from its instantaneous 
mean, where the barred quantities are the vortex positions av- 
eraged over z. To quantify the extent of vortex bending at a 
given time, we take the Fourier transform with respect to z 
(denoted by Fz) oiu{z^t), and calculate the power spectrum 



(10) 



Note that this measure of the power spectrum is insensitive to 
varicose waves |[T2|] and is thus entirely due to vortex bending. 

Fig. [3] shows the power spectra vs Kelvon wavenumber for 
a range of geometries at t = 0.70 s. We rescale the k axis 
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FIG. 5. Vortex decay rate as a function of relative temperature, for 
a range of geometries. Here the critical temperature is adjusted to 
account for interactions and finite size effects (Eq. (ED). The highest 
decay rate (for A = 20) is at T/T^^^^,,^ = 0.98, giving 1/i = 
3.9 s-\ 



by k = tt/W, which is the wavenumber corresponding to 
the largest visible wavelength within our choice of window 
length. In all geometries we see the power spectra take a 
thermalized form, with maximum occupation of the long- 
wavelength bending modes (low values of k). The thermal- 
ization of the vortex line from the straight initial conditions 
is most directly monitored using the Kelvon power spectrum. 
We find that typically within ~ 100 ms the Kelvon spectrum 
reaches a stable state. This is much shorter than the vortex 
decay time (see Sec. IIIIBI) . so we can rule out the influence 
of any initial state artifacts on our predictions for the vortex 
decay process. 



B. Effect of geometry: Kelvon power and vortex decay 

The effect of geometry on the importance of the Kelvon 
modes is clearly visible comparing Figs. [3](a)-(c): The total 
power in the spectrum decreases significantly as the system 
becomes more oblate. It is useful to calculate this total power 
in the bending modes, Ptot = / Pu{kj t)dk, which by Parse- 
val's theorem is related to the line-averaged deviation of the 
vortex from its mean position [|27|] . Results for the total power 
are shown in Fig.|4](a). It decreases exponentially as the ge- 
ometry becomes more two-dimensional, and we see that the 
total power is independent of oblateness for the most highly 
oblate condensates. We see that vortex bending is significantly 
reduced in the highly oblate systems, being almost entirely ab- 
sent for 8 < A. 

Fig. |4l;b) shows the effect of condensate geometry on the 
vortex lifetime, for a range of system temperatures. The life- 
time is quantified here in terms of the mean first exit time 1I22I1 . 
denoted by i. This is calculated as the ensemble average of the 



time it takes for the vortex to become indistinguishable from 
thermal fluctuations at the BEC edge (in the plane z = 0). The 
vortex decay rate, 1/f, decreases with increasing A, following 
the same trend as the total bending mode power [Fig. IHa)]. 
However, rather than decreasing to zero, the decay rate ap- 
proaches a constant value of order ~ 0.1s~^, being almost 
independent of A for 8 < A < 20. This constant is deter- 
mined by purely 2D vortex dissipation processes which do not 
involve bending. Interestingly, for the highest temperatures 
considered, the decay rate starts to increase again at A = 20. 
This increase is not associated with an increase of bending 
mode power [Fig.|Ha)]. Aside from this anomaly, discussed 
below, we infer a critical oblateness of Ac ^ 8 for the onset 
of 2DV dynamics in our system. 

C. Effect of temperature: Anomalous decay in quasi-2D 
systems 

The results of the previous section show that the vortex en- 
ters a 2D dynamics regime at A ~ 8, associated with strong 
suppression of bending modes, and a diminished vortex decay 
rate. However, for the highest temperatures the decay rates ex- 
hibit a revival as A increases. We interpret this as arising from 
the system entering a phase fluctuating quasi-2D regime. 

First, to more accurately characterize the critical tempera- 
ture we include the downward shifts from due to interac- 
tion and finite-size effects i28ll . This gives a corrected critical 
temperature 

where 5TI^ = -0.73|7V-^/^T^^ is the first-order correc- 
tion due to finite-size effects, and ST'J'^ = -1.33f iV^/^T^o 
is the first order correction due to interactions, with uj = 
{2ujr + ojz)/^ and a = {h/muj)^^'^ . In Fig. [5] we show the 
decay rates versus Tc(corr) verifying that these systems at the 
highest temperatures remain below the 3D critical tempera- 
ture. However, for oblateness 12 ^ A we see anomalously 
rapid increase in the vortex decay rate for temperatures ex- 
ceeding ~ 0.95Tc(corr)- For these cases the system crosses 
over to being quasi-2D as the chemical potential becomes 
of the order of the energy to activate degrees of freedom in 
the tightly confined direction, i.e. fi ^ hujz + (where 
eo is the single particle ground state energy). In particular, 
the highest temperature results for A = 8, 12, 15, 20 have 
l^/hujz = 2.2, 1.69, 1.42, 1.19 (c.f. huJz + eo ^ 3huJz/2). In 
the quasi-2D regime the dynamics of vortices are not well- 
understood, but strong phase fluctuations will likely facilitate 
the rate of vortex damping as we observe in our results [|29ll . 
A systematic investigation of vortex decay in a system more 
deeply in the 2D regime is therefore needed, but is beyond the 
scope of our work here. 

IV. DIPOLE DECAY 

Having systematically studied the role of geometry and 
temperature in single vortex decay, we now consider two par- 
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FIG. 6. Density isosurfaces of the C-field showing the SPGPE evolution of a vortex dipole with initial vortex locations (x, y±) = (0, ±s/2), 
for separation s — S/xm. Red: the C-field isosurface. Blue: isosurface restricted to the region of the vortex cores. 



ticular cases of the decay of a vortex dipole at finite tempera- 
ture. 



A. Spherical geometry: decay via ring and loop formation 



We first consider dipole decay in a spherical trap. 
We use SPGPE simulation parameters (T,/i, ecut) = 
(52 nK, 37ficj, 71 to), which give a total atom number N = 

2 X 10^ and dimensionless damping rate ^ of /ksT = 

3 X 10~^. These parameters give a system with the same 
atom number and temperature as the oblate configuration con- 
sidered in the experiment of Ref. [fisll and in the next sec- 
tion, so the system serves as a benchmark for the effect of 
geometry on vortex dipole dynamics in this system. We cre- 
ate a finite-temperature vortex dipole state by imposing a vor- 
tex dipole phase pattern on a SPGPE equilibrium state via 
V^c(x) ?/^c(x)e*®'^''\ where 6(x) = arctan (^^) - 
arctan (^^^)- This creates an axially aligned vortex dipole 
with positive and negative vortices located at = 
(0, ±5/2), with s = 3/im. Figure [6] shows the decay sequence 
for the vortex dipole. The initially axially aligned vortices 
rapidly curve as the vortices precess, such that the vortex cores 
are normal to the surface. This long- wavelength bending leads 
to formation of a vortex ring, which then subsequently disso- 
ciates into two vortices at t = 0.5s. As the vortices closely 
approach again at t = 0.62s, both loops and rings form, re- 
ducing the total vortex length, and leading to total decay in 
less than Is. We thus observe rapid vortex dipole decay via 
intermediate vortex ring and loop formation. 




FIG. 7. Average number of vortices as a function of time during vor- 
tex dipole decay for numerical simulations and experimental data. 
The dashed lines show SPGPE results for a range of initial vortex 
dipole separations. Circles show experimental data of Ref. |15]. The 
blue diamond at t = s shows the average number of vortices ob- 
served in the experimental initial state due to spontaneous vortex for- 
mation Q. 



B. 2DV regime: comparison of experiment and theory 

In this section we use ([5]) to model the dynamics of a 
charge-one vortex dipole in a highly oblate finite-temperature 
Bose-Einstein condensate, as in the experiment by Neely et 
al [15]. The system consists of ~ 2 x 10^ atoms at tem- 
perature ~ 52nK, held in a harmonic trap with frequencies 
{ujr.^z) = 27r X (8,90)Hz. The trap aspect ratio A = 11.25 
exceeds Ac ~ 8 identified in Sec. IIIIBI for closely related 
parameters. We can therefore expect that the experimental 
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FIG. 8. Column densities (top) and phase slices for the C-field showing the SPGPE evolution of a vortex dipole with initial vortex locations 
(a:, y±) — (0, =bs/2), for separation s = 3 \mi. The series spans two orbits of the dipole, where a single orbital period is ^ 1.2s, culminating 
in mutual vortex annihilation towards the completion of the second orbit. The peak density is renormalized to unity at each time (upper 
colorbar), and the phase is represented on the interval [— tt, tt) (lower colorbar). 



regime is also approximately that of 2DV. Experimentally, a 
single vortex dipole is created by pushing the BEC around an 
past an obstacle above a critical velocity. Ideally, we would 
like to compare the distance s between the two vortices upon 
the creation of the dipole in both experiment and simulations, 
and to relate s to the vortex dipole lifetime in each case. How- 
ever, the experimental spatial resolution was ~ 6/im, which 
is also the approximate vortex separation found immediately 
after dipole creation in Gross-Pitaevskii simulations of the 
expriment |[l5|] . Since the vortices are not individually re- 
solvable immediately after creation, the experimental data can 
only place an upper bound on s. We are thus unable to carry 
out quantitative comparison at this point. In this work we in- 
stead map out the dependence of dipole lifetime on s by com- 
puting the decay dynamics for a range of initial vortex sepa- 
rations. 

As in the spherical geometry of Sec. IIV Al we use SPGPE 
simulation parameters (T,/i, ecut) = (52 nK, 37ficD, Tlfio)), 
which again gives a total atom number = 2 x 10^ and di- 
mensionless damping rate (|9]) of fi'^/kBT = 3 x 10~^. As 
in the spherical case, we phase imprint a vortex dipole into 
the BEC. We choose the initial separation between vortices as 
s = 1,2, 2.5, 3, and 3.5 /im. Note that the condensate radius 
is ~ 50/im, so these separations are small relative to the sys- 
tem size, but larger than the healing length ^ ^ 0.2/im. Qual- 
itatively, our simulations differ from the experiment in several 
ways: in our simulations we have two vortices initially, a high 
degree of cylindrical and mirror symmetry, and a well defined 
total atom number. In the experiment there is occasionally a 
vortex present from the BEC formation process [|2Tll , an irre- 
ducible (albeit very small) cylindrical asymmetry of the har- 
monic potential, a mirror asymmetry of the stirring potential. 



and uncertainty in the total atom number. We do not attempt 
to model all of this complexity here, confining our study to the 
basic decay process of internal annihilation in the symmetric 
system, with fixed total atom number ||3Q|i . 

Figure [7] shows the mean vortex number over time from the 
experimental data |[T5|] . compared with the results of our sim- 
ulations. In our simulations we always find that the vortices 
mutually annihilate each other in the center of the BEC, rather 
than damping at the boundary. The vortices move rapidly 
when closely approaching each other, nearing the speed of 
sound c = 0.2cm/s. Consequently, most of the orbital time 
is spent at the outer edge of the BEC, and the time interval 
in which annihilation may occur is a small fraction of the or- 
bital period (~ 10%). Thus the vortex number has an abrupt 
time dependence for a given initial separation, with character- 
istic timescale given by the orbital period of the vortex dipole, 
found numerically and experimentally to be ~ 1.2s and inde- 
pendent of s over a wide range [15.1 . In general, the s value 
selects a given orbit of annihilation, within which there is very 
small variation of lifetime. 

In detail, we find that for 5 = (1, 2, 2.5, 3, 3.5)/imthe mean 
lifetime is r = (0.01, 0.86, 2.04, 2.07, 3.30)s, corresponding 
to A'o = (O7 2, 2, 3) dipole orbits before annihilation. For 
example, with s = 1/im, the dipole self annihilates immedi- 
ately after creation. The separations s = 2.5/im and s = 3fim 
lead to almost identical lifetimes of about 2s, corresponding to 
the significant drop in vortex number seen experimentally for 
2s ^ t ^ 2.5s which we interpret as an experimental indica- 
tion of the dipole lifetime. 

In Figure [8] we plot a time series of simulation data for 
s = 3/im showing internal annihilation at the end of the sec- 
ond orbit. Interestingly, we observe a tilt in the axis of dipole 
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precession caused by thermal fluctuations, which is also ob- 
served in experimental absorption images. 

V. CONCLUSIONS AND OUTLOOK 

Here we summarize our results on single vortex and vortex 
dipole decay. 

Sub-critical Regime. — For spherical systems at T = 
0.78T^^, a lone vortex is freely deformed from its initial 
straight-line configuration by the activation of vortex bend- 
ing modes. This effect diminishes as the condensate is flat- 
tened, and bending modes become energetically inaccessible. 
The total bending mode power decreases exponentially as the 
condensate geometry becomes increasingly two-dimensional. 
For geometries with oblateness A > Ac(^ 8), the total power 
asymptotes to zero, indicating bending is strongly inhibited. 
The effective vortex decay rate also decreases as A increases, 
approaching a constant two-dimensional value for Ac < A. 
Ac can therefore be identified as the critical oblateness re- 
quired for the system to enter the 2D vortex regime. 

Critical Regime. — The foregoing comments hold for tem- 
peratures outside an observed critical regime. When critical 
fluctuations associated with the phase transition become sig- 
nificant, we find that the vortex decay rate is highly sensi- 
tive to temperature and geometry. In particular, we see an 



anomalously large increase in the decay rate for geometries 
with Ac < A, and 0.95Tc(corr) ^ T that we attribute to phase- 
fluctuations associated with the system entering the quasi-2D 
regime. 

Vortex Dipole Decay. — We have modeled the experiment of 
Neely et al which is in the 2DV regime (A = 11.25), and 
found that the approximate time for decay of a vortex dipole 
in the experiment (~ 2.25s) is consistent with numerical sim- 
ulations with an initial vortex separation of s 2.5 — 3/im. 
In the simulations with these initial conditions mutual vortex 
annihilation occurs after two complete orbits, at t ~ 2s. In 
this work we have used s as the fitting parameter; improved 
experimental resolution would be required for a quantitative 
comparison without any fitting parameters. In similar compu- 
tations of spherical (A = 1) BECs, the vortices decay much 
more quickly to rings and other bent geometries that would 
not be readily visible after only ~ 100ms using the experi- 
mental methods of Ref . ifTsIl . We thus infer from the exper- 
iments of Ref. fisl] and the calculations presented here that 
the experimental system consisting of a BEC with oblateness 
A ^ 11.25 is within the 2D vortex dynamics regime, consis- 
tent with our calculations suggesting a critical oblateness of 
Ac ^8. 
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